home *** CD-ROM | disk | FTP | other *** search
/ Fritz: All Fritz / All Fritz.zip / All Fritz / FILES / EDUCNOMY / ASTROSET.LZH / SUNSHI.BAS < prev    next >
BASIC Source File  |  1985-12-20  |  896b  |  40 lines
  1.  
  2.  
  3. 10 REM    SUNSHINE
  4. 12 REM
  5. 14 P=3.14159265: R=P/180
  6. 16 INPUT "LATITUDE (DEGREES)";F
  7. 18 F=F*R
  8. 20 S=SIN(F): C=COS(F)
  9. 22 INPUT "START DAY IN 1986 ";J1
  10. 24 INPUT "END DAY IN 1986   ";J2
  11. 26 IF J2<J1 THEN 24
  12. 28 Z=0
  13. 30 FOR D=J1 TO J2
  14. 32 M=357.1663+0.9856*D
  15. 34 M=M*R
  16. 36 L=279.8603+0.985647*D
  17. 38 L=L+1.91533*SIN(M)
  18. 40 L=L+0.02001*SIN(2*M)
  19. 42 L=L+0.00029*SIN(3*M)
  20. 44 L=L*R
  21. 46 G=0.3978062*SIN(L)
  22. 48 K=SQR(1-G*G)
  23. 50 U=(-0.01454-S*G)/(K*C)
  24. 52 IF ABS(U)<1 THEN 58
  25. 54 IF U>1 THEN 70
  26. 56 Z=Z+1440: GOTO 70
  27. 58 F=0.9043+0.03023*COS(M)
  28. 60 F=F/(1-0.15825*SIN(L)*SIN(L))
  29. 62 A=ATN(SQR(1-U*U)/U)
  30. 64  A=A/R
  31. 66 IF U<0 THEN A=A+180
  32. 68 Z=Z+2880*A/(360.985647-F)
  33. 70 NEXT D
  34. 72 PRINT
  35. 74 PRINT "SUNSHINE = "
  36. 76 PRINT INT(Z+0.5);" MINUTES"
  37. 78 PRINT INT((Z+30)/60);" HOURS"
  38. 80 END
  39.  
  40.